BaFe 2 Se 2 as an Iron-Based Mott Insulator with Antiferromagnetic Order 
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A iron-based compound with a quasi-two-dimensional array of FeSesO tetrahedra and an or- 
thorhombic structure, namely BaFe2Se20, has been successfully fabricated. Experimental results 
show that this compound is an insulator and has an antiferromagnetic (AF) transition at 240 K. 
Band structure calculation reveals the narrowing of Fe 3d bands near the Fermi energy, which leads 
to the localization of magnetism and the Mott insulating behavior. The large distances between the 
Fe atoms perhaps are responsible for the characters. Linear response calculation further indicates 
a strong in-plane AF exchange J, this can account for the enhanced magnetic susceptibility (which 
has a maximum at about 450 K) above the Neel temperature. 
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A. Introduction 

Discovery of the copper-based superconductors^ with 
high transition temperatures (T c ) has attracted great at- 
tention and triggered extensive research in the fields of 
condensed matter physics and material sciences. High-T c 
superconductivity in the cuprates occurs with either hole 
or electron doping in the antiferromagnetic (AF) Mott in- 
sulating parent compounds. The Mottness clearly plays 
an important role here^. In the past years, scientists have 
been trying to explore new high-T c superconductors in 
non-copper-based compounds. Until four years ago, the 
discovery of superconductivity at 26 K in oxy-arsenide 
LaFeAsOi-^F^ broke the uniqueness of cuprates in giv- 
ing rise to high-T c superconductivity. This is actually a 
surprise since Fe atoms are normally strongly magnetic in 
materials and usually inhibit superconductivity. Among 
the iron-based superconductors with several different 
structures 3 - - — , it is widely perceived that the common 
FeX (X — pnictide or chalcogenide) layers are respon- 
sible for the superconductivity. Meanwhile, for many 
FeAs-based superconductors, their parent phases possess 
an AF order which seems to compete with the supercon- 
ductivity. To unravel the mystery of the origin of mag- 
netism and its relation with superconductivity is the most 
important issue. Now there exist two major pictures, 
namely the electron pairing is established through a weak 
coupling by exchanging the AF spin fluct uat ion s 1 2 1 1 3 , or 
the pairing is furnished by the local strong pairin g 14 ' 15 . 
Related to this debate, it has been argued whether the 
magnetism in the iron-based superconductors is local- 
ized or itinerant in nature. It thus becomes very im- 
portant to find a iron-based compound with a formally 
known knowledge of localized or itinerant magnetism. 
Recently, exploration along this line was successful in 
finding the materials of J Re 2 2 Fe20(Se/S)^^ (Re = 
rare earth element), La2Co2Se203^, L^OaM^Se^ 
and Na2Fe2Se202£. These compounds all contain an ex- 



panded similar square-planar M lattice in their M20Se2 
(M = Fe, Co and Mn) layers. Here we report the fabri- 
cation of an iron-based compound BaFe2Se20. From the 
structure point of view, it looks like the existing materials 
i?e202M 2 0(Se/S)2 and Na2Fe2Se20. However, a closer 
scrutiny finds the clear difference between them, since the 
compound BaFe2Se2 contains a honeycomb-like Fe sub- 
lattice instead of the square-planar Fe lattice. The com- 
pound has an orthorhombic structure and is revealed to 
be a Mott insulator with a localized antiferromagnetism. 



B. Experimental Details 



The single-crystalline sample was synthesized with the 
starting materials BaO and FeSe in the ratio of 1:2 by 
using a flux method, similar to a previous synthesis of 
K x Fe2- y Se-2r^. The resultant crystals were black shin- 
ing platelets with a typical dimension of 1.5 x 1.5 x 0.5 
mm 3 . To determine the compositional ratio of atoms in 
the compound, scanning electron microscope and energy- 
dispersive x-ray spectroscopy (EDX) measurements were 
carried out on the crystals. Powder X-ray diffraction 
(XRD) pattern was measured on the smashed and pow- 
dered crystals with 8 — 20 scan at room temperature. 
The XRD data were analyzed by a Rietveld method with 
the software GSAS 22 . Resistivity measurement was done 
using a four-probe technique on the physical property 
measurement system (PPMS). Dc magnetic susceptibil- 
ity was measured on a superconducting quantum inter- 
ference device (SQUID). Specific heat was also measured 
using a relaxation technique on the PPMS. 
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C. Experimental Data and Discussion 

1. Crystal Structure 

The energy-dispersive x-ray spectroscopy (EDX) of 
the crystals is shown in Fig. 1(a). One can see that 
the crystals are composed of Ba, Fe, Se and O ele- 
ments. The ratio of atoms is determined as Ba:Fe:Se = 
20.61:39.94:39.45 (close to 1:2:2) while the amounts of O 
can not be quantified accurately. The EDX result offers 
us a basis for seeking suitable reference compound to do 
phase identification. We finalize our compound isostruc- 
tural to /3-BaFe 2 S 3 (formerly reported as Ba 2 Fe4Sg22i2i). 
/3-BaFc 2 S3 may be viewed as an average due to twin- 
ning of a polymorph of ev-BaFe 2 S3 which contains iron- 
deficient-FeAs-layers-like Fe 2 S3 layers^—. BaFe 2 Se 2 
can be considered as a derivative of /3-BaFe 2 S3 by replac- 
ing the S atoms at 2b sites with O atoms and the S atoms 
at 4e sites with Se atoms. The powder x-ray diffraction 
(XRD) pattern and Rietveld fit for the smashed and pow- 
dered crystals is shown in Fig. 1(b). The data is well 
fitted with an orthorhombic BaFe 2 Se 2 phase which has 
the space group of Pmmn. The unit cell is given in the 
inset of Fig. 1(b). The starting parameters for the fit 
are taken from /?-BaFe 2 S3^ and the GSAS program fi- 
nally finds the best fitting parameters. In Table I, the 
best fitting structure parameters are listed. As shown 
in Fig. 1(c), in BaFe 2 Se 2 0, the nearest neighbor atoms 
of each Fe atom, namely the three Se atoms and one O 
atom, surround the Fe atom and form a tetrahedra. The 
structure of BaFe 2 Se 2 can also be viewed as an infinite 
quasi- two-dimensional array of FeSe30 tetrahedra. Up- 
per and lower Fe ladders along the a direction and zig-zag 
Fe chains along the b direction forms uneven Fe layers. 
Obviously, there is not a strictly speaking square-planar 
Fe lattice in this compound. It is interesting to notice 
that the upper and lower Fe ladders form a non-coplanar 
honeycomb- like structure, as shown in Fig. 1(d). 



2. Transport properties 

Fig. 2 gives the resistivity measurement result for 
the BaFe 2 Se 2 single-crystalline sample. The room- 
temperature resistivity measured is as high as 100 f2 
cm which is comparable to those of La 2 2 Fe 2 OSe2^, 
Na 2 Fe 2 Se 2 OS2 and Ba 2 F 2 Fe 2 OSe2^. In the tempera- 
ture dependence of resistivity, an insulating behavior is 
observed, and the resistivity divergence in the low tem- 
perature region below 150 K can be fitted by the model 
of variable-range- hopping: p(T) — poexp(To/T)i . This 
result was also observed in many 3d transition metal com- 
pounds in which the correlation is assume to be strong. 

The temperature dependence of dc magnetic suscep- 
tibility for the BaFe 2 Se 2 single-crystalline sample is 
presented in Fig. 3(a). The dc magnetic susceptibility 
measurement was carried out under a magnetic field of 
20000 Oe in zero-field cooling (ZFC) process. We can see 
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FIG. 1: (Color online) (a) Energy-dispersive x-ray spec- 
troscopy (EDX) of the crystal. The inset shows a crystal with 
the typical dimension of 1.5 x 1.5 x 0.5 mm 3 , (b) Powder 
X-ray powder diffraction (XRD) pattern and Rietveld fit for 
the smashed and powdered crystals. The inset shows a unit 
cell of BaFe 2 Se 2 0. (c) Crystal Structure of BaFe 2 Se 2 0. (d) 
Honeycomb- like Fe sub-lattice in BaFe 2 Se 2 0, the gray balls 
and the red balls all indicate Fe atoms, with the gray balls 
forming the upper ladder and the red balls forming the lower 
ladder, respectively. 



TABLE I: Crystallographic data of BaFe 2 Se 2 Q. 



Space group 
a (A), 6 (A), c (A) 

V(A 3 ) 
z 

R p , wR p 
Atomic parameters: 
Ba 

Fe 

Se 

O 

Bonding lengths (A): 
Ba-O 
Fe-0 
Fe-Se 
Fe-Se 
Fe-Fe 
Fe-Fe 

Bonding angles (degree) 
Se-Fe-Se 
Se-Fe-Se 
Se-Fe-0 
Se-Fe-0 



Pmmn 

4.14028(5), 9.86716(14), 6.73381(8) 
275.095 
2 

2.41%, 3.12% 

2a (i |, z) 
z=0.50850(18) 
4e (i, y 7 z) 
y=0.91347(22),' z=0.87920(27) 

4e (i, y 7 z) 
y=0.95840(15),' z=0.24251(20) 
2b (|, |, z) 
2=0.7502(13) 

2.706(6) x 2 

1.832(4) xl 
2.5604(15) x2 
2.4863(23) xl 
3.1382(13)x2 
3.2260(12) xl 

103.11(6)x2 
107.90(9) xl 
106.44(14) x 2 
128.57(29) xl 
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FIG. 2: (Color online) Temperature dependence of resistivity 
for the BaFe2Se20 single-crystalline sample. The inset shows 
that hip is directly proportional to J 1-1 / 4 below 150 K. 



a kink associated with the AF transition at about 240 
K in the x(^) curve. This transition is indicated more 
clearly by the peak in the d(xT)/dT versus T curve, 
as shown in the inset of Fig. 3(a). Similar AF transi- 
tion is observed in the parent compounds of cuprates and 
some iron-based superconductors 2 ^— . Above 240 K, the 
magnetic susceptibility exhibits a continuing enhance- 
ment, which is observed in many two-dimensional mag- 
netic systems including the parent compounds of iron- 
based superconductor o ' 31 and is explained as due to the 
short-range correlation of the local moments 3 . 2 . The con- 
tinuing enhancement leads to a broad maximum around 
T m ax = 450 K. The ratio of Tpj/T max is usually used to 
estimate the extent of the low-dimensional magnetic cor- 
relations. As Tjv /T max in our BaFe2Se20 is 0.53, it can 
be concluded that BaFe2Se20 is a quasi- two-dimensinal 
magnetic system. 

In order to further understand the magnetic transition, 
we performed a specific heat measurement, as shown in 
Fig. 3(b). Provided the negligible contribution by the 
magnons at very low temperatures (T <C Tjy), the spe- 
cific heat can be separated into electronic and phonon 
parts: C(T) = jT + /3T 3 . We find that the low tem- 
perature data show no electronic contribution, consistent 
with the insulating ground state. From the slope (/?) of 
the straight line in the C/T vs T 2 plot (as shown in the 
inset of Fig. 3(b)), we obtain the Debye temperature ®d 
= 224.5 K using the formula 9^ = [12tt 4 NR/ (5 0)]* , 
where N is the number of atoms in the chemical formula 
and R refers to the ideal gas constant 8.314 J mol -1 K _1 . 
Since the Debye model fails to reproduce the C(T) data 
at high temperatures, we employ a polynomial to fit the 
phonon contribution in the temperature range from 50 to 
300 K. Assuming that the total specific heat consists of 
phonon and magnetic components, the magnetic contri- 
bution can be obtained simply by the subtraction of the 
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FIG. 3: (Color online) (a) Temperature dependence of dc 
magnetic susceptibility for the BaFe2Se2 single-crystalline 
sample in zero-field-cooling (ZFC) process at a magnetic field 
of H = 20000 Oe. The inset shows the d(xT)/dT versus T 
curve. We can see a peak associated with an AF transition at 
240 K. (b) Temperature dependence of specific heat for the 
BaFe2Se20 single-crystalline sample. The green solid curve 
represents the phonon contribution, fitted by a polynomial. 
After deducting the phonon contribution, the specific heat 
contributed by magnetic shows a peak at 240 K. Part of the 
magnetic entropy has been dissipated above 240 K. 



phonon contribution. A heat-capacity peak appears at 
240 K which is exactly the temperature where the mag- 
netic susceptibility exhibits a sharp step (as illustrated 
by the derivative d(%T) /dT) . The magnetic entropy can 
be calculated using the integral S mag = J C mag /T dT. 
The result gives S mag = 4.89 J mol" 1 K _1 . Part of the 
magnetic entropy has been dissipated above 240 K. Al- 
though the AF occurs at about 240 K, as reflected from 
the d(%T)/dT data, we find a significant contribution of 
entropy due to the strong AF fluctuation above 240 K. 



D. Band Structure Calculations 

To have a comprehensive understanding on the physi- 
cal properties of BaFe2Se20, we also perform the density- 
functional theory calculations based on the local density 
approximation (LDA) to density functional theory with 
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FIG. 4: Band structure from LDA calculation, and LDA+ U 
calculation with U = 6.0 eV. 



the full-potential, all-electron, linear-muffin-tin-orbital 
method 33 . To see the basic features of the electronic 
structure, we first do nonmagnetic LDA calculation, and 
show the band structure in Fig. 4a. Although with 
strong hybiridization, the Fe, Se and O bands are sep- 
arated as shown in Fig. 4a. The six bands located 
from -6.6 to -5.0 eV basically come from O 2p states, 
the twelve bands located from -5.0 to -3.4 eV are mainly 
contribuated by Se 4p states, while the twenty Fe 3d 
bands distribute from -1.4 to 1.2 eV. The non-coplanar 
honeycomb- like structure significantly narrows the band, 
consequently, the bandwidth of Fe 3d, Se Ap and O 2p 
bands are considerablely narrower than these of Fe-based 
superconductor and La2C>2Fe20Se2 16 . 

We then perform the spin polarized calculation. We 
first evaluate the interlayer magnetic exchange inter- 
action based on our linear response scheme 3 ^. Our 
numerical results show that similar with Fe-based 
superconductor-^ 5 ., the exchange process happens almost 
completely in the Fe-Se-0 layer, and the inter-layer ex- 
change interaction is negligible. This result is consistent 
with the observed two-dimensional magnetic behavior. 
To search the ground state magnetic ordering configu- 
ration, we consider four spin alignments within layer as 
following: FM, AFM-1 (Fel and Fe2, Fe3 and Fe4 couple 
ferromagnetically while Fe2 and Fe3 have opposite spin 
orientation), AFM-2 (Fel and Fe2, Fe3 and Fe4 couple 
antiferromagnetically while Fe2 and Fe3 have the same 
spin orientation), AFM-3 (Fe atom couples with all of 
three nearest neighbor Fe antiferromagnetically). The 
three kinds of AFM spin alignments are plotted in Fig. 
5. The LSD A calculation predicts that the AFM-3 has 
the lowest total energy, and the magnetic moment mainly 
locates at Fe site. It is interesting to notice that the cal- 
culated magnetic moment (about 3.0 (j,b) is not sensitive 
to the magnetic configuration. While for iron arsenides, 
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FIG. 5: Three kinds of possible AFM spin alignments within 
layer. AFM-3 has the lowest total energy. 



the magnetism is quite itinerant and the magnetic mo- 
ment strongly dependent on the magnetic configuration. 
Although can explain the AFM behavior, LSDA calcula- 
tion fails to reproduce the observed insulating behavior. 

It is well known that the Coulomb interaction for 3d 
orbital is significant, and estimates for the values of U 
have been recently done yielding U rs 5.0 eV in Fe-based 
superconductors 3 -^. Although the accurate value of U 
is not known for BaFe2Se20, with narrower bands we 
generally expect it is larger than that of Fe-based su- 
perconductor. We therefore utilize LDA+ U scheme 3 -^, 
which is adequate for the magnetically ordered insulating 
ground states 3 -^, and set U = 6.0 eV. Same with LSDA 
calculation, LDA+ U calculation also predicts AFM-3 as 
the ground state, and the obtained magnetic moment is 
not depending on the magnetic configuration. As shown 
in Fig. 4b, electronic correlation significantly changes the 
band structure and results in a considerable effective gap 
(about 1.86 eV). Our linear response calculation 35 shows 
that the magnetic exchange interaction is short range 
and the magnetic coupling further than the second near- 
est neighbor is almost equal to zero. The magnetic ex- 
change interactions between all of three nearest neigh- 
bor Fe-Fe are AFM sign, which is consistent with the 
total energy comparison. The strength of exchange cou- 
pling {J i = 96 meV and Ji = 31 meV) is similar with 
that of LaFeAsO 35 ' 39 , which may account for the fact 
of similar magnetic susceptibility above the Neel tem- 
perature. We also vary parameter U between 2 and 8 
eV for Fe 3d electrons to see what effects the on-site 
Coulomb repulsion would bring to the electronic struc- 
ture of BaFe2Se20. Our numerical calculation confirm 
that U does not change the qualitative results: AFM- 
3 is the ground state, the magnetism is local in nature, 
and the Coulomb interaction is essential for opening the 
band gap. The real magnetic structure will be deter- 
mined by neutron diffraction method, just like the case 
of AeaOaFezSeziLaa. 



E. Concluding Remarks 

In summary, a iron-based compound BaFe2Se20 was 
synthesized successfully. It has an orthorhombic struc- 
ture with the space group of Pmmn. Strictly speaking, 
there is not a square-planar Fe lattice in it, instead, it 
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contains a honeycomb-like Fe sub-lattice. The model of 
variable-range-hopping is used to explain the resistivity 
divergence in the low temperature region. An AF or- 
dering transition at 240 K is observed both from the dc 
magnetic susceptibility and specific heat. Band structure 
calculation reveals the narrowing of Fe 3d bands near 
the Fermi energy, which leads to the localization of mag- 
netism and the Mott insulating behavior. The distances 
between the Fe atoms (3.1382(13) A and 3.2260(12) A) 
perhaps are responsible for the characters, since they are 
so large that the overlap between the 3d-orbitals becomes 
very small. Linear response calculation further indicates 
a strong in-plane AF exchange J, this can account for 
the enhanced magnetic susceptibility (which has a maxi- 
mum at about 450 K) above the Neel temperature. Con- 
sidering the fact that the Mottness and AF order are 



commonly occurring in cuprates and the iron pnictides, 
it would be interesting to know whether it is possible 
to induce superconductivity by either chemical doping 
or applying high pressure in our present system. This 
investigation is underway. 
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